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Consistent distances from Baade-Wesselink analyses of 
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ABSTRACT 

By using the same algorithm in the Baade-Wesselink analyses of Galactic RR Lyrae 
and Cepheid variables, it is shown that, within 0.03 mag ler statistical error, they yield 
the same distance modulus for the Large Magellanic Cloud. By fixing the zero point 
of the color-temperature calibration to those of the current infrared flux methods and 
using updated period-luminosity-color relations, we get an average value of 18.55 for 
the true distance modulus of the LMC. 
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1 INTRODUCTION 

A large number of papers have been published lately on the 
determination of the distance to the Large Magellanic Cloud 
(LMC). This increased interest is attributed to the extreme 
values published after the analyses of the HIPPARCOS data 
on Cepheids (Feast & Catchpole 1997) and to the applica- 
tion of the red clump stars of the LMC through the same 
type of Galactic stars calibrated by HIPPARCOS (Stanek et 
al. 2000) . Along this line, there seems to exist a general be- 
lief that Cepheids yield usually longer distances (distance 
moduli of about 18.5 mag and above), whereas RR Lyrae 
stars give shorter ones (~ 18.3 mag). This latter belief is 
based on statistical parallaxes (Popowski & Gould 1998), 
Baade-Wesselink (BW) analyses (see Clementini et al. 1995, 
and references therein, hereafter C95) and the application 
of the HIPPARCOS data (Fernley et al. 1998; see however 
Groenewegen & Salaris 1999). It is important to note that 
distances derived from the statistical parallax method em- 
ploy rather small sample of stars (« 200 objects), whereas 
the HIPPARCOS parallaxes for RR Lyrae stars are very in- 
accurate. Actually, more precise HST parallaxes yield dis- 
tances in good agreement with the "canonical" LMC value 
of 18.5 mag derived from both prototype variables RR Lyr 
and S Cep (Benedict et al. 2002a,b). It is also important 
to mention that all results based on the pulsation analyses 
of double-mode variables (RR Lyrae, Cepheids, both funda- 
mental & first overtone and first overtone & second overtone 
ones) give basically the canonical value for the LMC distance 
(see Kovacs 2002 and references therein). 

In the light of the above, it is important to address 
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the question of consistency between the BW analyses of 
Cepheids and RR Lyrae stars. Previous works have been per- 
formed by different groups of researchers, often using differ- 
ent BW parameters (e.g., p factors, the ratio of the observed 
and surface radial velocities) and effective temperature rela- 
tions (surface brightness for Cepheids - e.g., Laney & Stobie 
1994, stellar atmosphere models for RR Lyrae stars - e.g., 
Liu & Janes 1990). In addition, there is a role of the some- 
what delicate nature of the numerical method employed in 
the actual implementation of the BW analysis (e.g., compu- 
tation of the surface acceleration). Therefore, we decided to 
visit the problem by using the same BW algorithm for both 
groups of stars. In order to compute the LMC distance mod- 
ulus, the derived absolute magnitudes will be used in con- 
junction with currently published Period-Luminosity-Color 
(PLC) relations (Udalski et al. 1999; Kovacs & Walker 2001, 
hereafter U99 and KW01, respectively). 



2 IMPLEMENTATION OF THE BW METHOD 

Here we basically follow the standard methodology of the 
BW analysis as usually implemented in RR Lyrae studies 
(e.g., Liu & Janes 1990). In brief, we use the observed V 
light, V — K color and Kad radial velocity curves to compute 
the luminosity (L), effective temperature (T e a) and radius 
(R) variations. For any assumed value of the equilibrium 
test radius Ro we compute the following quantities: 



Temporal radius: 

- R{t) = R + AR(t), 

- AR(t) = - P J*[v iad ( 
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where Pq denotes the period, p is the velocity projection 
factor, that we assume to be constant and set equal to 1.35. 
(We note that Gieren, Fouque & Gomez 1997 used a period- 
dependent p factor. Their values span the range of 1.35- 
1.37 in the period regime of our sample. The effect of using 
lower p is minute, because an increase of 0.02 in p yields an 
absolute magnitude that is brighter only by 0.02 mag.) 
• Temporal T e g : 

- logT off = f[(V - K) , [Fe/H],log fl ], 

where /[...] denotes the interpolation polynomial of the stel- 
lar atmosphere models of Castelli, Gratton & Kurucz (1997) 
with the zero point determined by the scale of Blackwell 
& Lynas-Gray, (1994) (hereafter BLG94). The models have 
been computed without convective overshooting. The mi- 
croturbulent velocity has been fixed at 2 kms -1 . The distri- 
bution of heavy elements follows that of the Sun. In order 
to find the temperature zero point offset between the model 
values of Castelli et al. (1997) and the ones given by the 
infrared flux method of BLG94, the compilation of the ob- 
servational data on non-variable stars by C95 is used. We 
get log T cff (model)- log T cff (BLG94)= 0.0036. 

It is interesting to compare this temperature scale with 
the one implied by the surface brightness (SB) method 
of Fouque & Gieren (1997). Rewriting their formula given 
for the apparent angular diameter, we get logT c ff = 
-0.131(U - K) - 0.1BC + 3.9474. By using the above 
mentioned compilation of C95, we get an average difference 
of logT eff (SB)-logT eff (model& BLG94)= -0.0035. By per- 
forming the same test on the BW samples (see Table 1), the 
differences are positive, but less than ~ 0.002 for Cepheids. 
For RR Lyrae stars, the differences are even smaller, with 
changing signs. The larger negative difference on the sample 
of C95 is due in part to the larger weight of the higher 
temperature stars. In conclusion, we may state that the 
present temperature scale in the temperature range of in- 
terest is very similar to the ones used in earlier BW studies 
of Cepheids. 

For RR Lyrae stars, the iron abundance [Fe/H] is com- 
puted from the Fourier decompositions of the V light curves 
as given by Jurcsik & Kovacs (1996). For Cepheids, we set 
[Fe/H]=0, which seems to be a reasonable first approxi- 
mation (see Luck et al. 2003 for the metallicity distribu- 
tion of Cepheids). The temporal gravity g is computed by 
g = GM/R 2 — p dVrzd/dt, where G is the gravitational con- 
stant, M is the stellar mass, which is fixed at 0.7Mq and 
5.0Mq for RR Lyrae stars and Cepheids, respectively. Since 
the contribution of log g to log T c g is relatively small when 
the V—K color index is used, the exact value of M is not cru- 
cial. In principle, we can compute a consistent value for M 
by using the pulsation equation (e.g., van Albada & Baker 
1973). However, both the weak dependence of logT of f on 
logg, and the difference between the equilibrium and static 
stellar parameters (see Bono, Caputo & Stellingwerf 1995) 
make this procedure less straightforward and therefore, is 
left for the task of future studies. The dereddened color in- 
dex (V — K)o is computed from the observed one by us- 
ing Eb-v reddening values either as given in the Galactic 
Cepheid Database of the David Dunlap Observatory 1 or by 

1 http://ddo.astro.utoronto.ca/cepheids.html 



the corresponding formula of KW01 for RR Lyrae stars. Be- 
cause of the large reddenings of the Galactic Cepheids, we 
use the formulae for selective absorption ratios as given by 
the BW studies of Cepheids (see, e.g., Gieren, Fouque & 
Gomez 1998): 

- R v = 3.07 + 0.28(B - V)o + 0.04E B -v, 

- Ri = 1.82 + 0.205(B - V) + 0.022E B -v, 

- Rk = 0.279, 

- Av = RvEb-v, Ai = RiEb-v, Ak = RkEb-v , 

where the / and K colors are in the Kron-Cousins and stan- 
dard Johnson systems, respectively. For compatibility, these 
reddening relations are used also for RR Lyrae stars. 

• Temporal absolute V magnitude: 

- M v = -2.51ogL + M b ©-SC(logTeff,[Fe/H],log 5 ), 

- log L = 4 log T cH /T® + 2 log R/Rq, 

where Rq and T c ® are the solar values. The latter one and 

the solar bolometric magnitude M® ol are set equal to 5778 K 
and 4.75 mag, respectively. The temporal bolometric correc- 
tion BC is computed from stellar atmosphere models in the 
same way as logT c ff above. The zero point of BC is set to 
a value which yields BCq = —0.07. 

• Point-by-point and average distance moduli: 

- DM(i) = V(i) - R v E B -v(i) - M v (i), 

- (DM}=N- 1 J2^ 1 DM(i), 

- o- 2 DM = (N- l)- 1 Eli[DM(i) - (DM)] 2 , 

where N is the number of data points on the light curve. 

The best fit is obtained at the test radius Ro which 
yields the smallest a dm- Usually the fits are fairly accu- 
rate, with a dm = 0.02-0.01. It is also worth mentioning the 
following technical details: 

• When straightforward Fourier fitting does not work, 
smooth light and velocity curves are obtained through 
Fourier decompositions of the high order polynomials fitted 
by least squares to the observed data points of the folded 
light curves. 

• Tables of stellar atmosphere models are employed 
through quadratic interpolation. 

• All pulsation phases are used, including also the possi- 
bly shock-disturbed quick expansion phase. 

• In spite of the careful treatment of the numerical deriva- 
tion of the radial velocity curve, the total gravitational ac- 
celeration occasionally reaches values which lie outside the 
tabulated values. In these brief moments of pulsation the 
closest extreme values as given in the tables are used. The 
same method is employed for filtering out the non-tabulated 
values of T e g. We note that these out-of-range errors occur 
almost always out of the minimum of a dm, where the as- 
sumed unphysical stellar radius further enhances the numer- 
ical problems related to logg. 



3 PLC RELATIONS 

Once the optimum equilibrium stellar parameters are com- 
puted by the above BW technique, we can evaluate the 
zero points of the corresponding PLC relations. Both for 
Cepheids and for RR Lyrae stars the PLC relations have 
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been derived from stellar populations different from the 
Galactic BW samples. In the case of RR Lyrae stars there 
is a reasonable large overlap between the metallicities of 
the present BW sample and the cluster variables used by 
KW01 to derive the empirical formula. Therefore, we may 
assume that the formula may also be valid for most of the 
field variables used in the BW analysis. There is also a dis- 
pute on the possible metal dependence of the zero point of 
the Cepheid PLC relation. However, we think that there 
are more theoretical and empirical results which indicate 
the non-existence of such a dependence (Alibert et al. 1999; 
Bono et al. 2002; Turner & Burke 2002). 

In the case of Cepheids we need to evaluate Cvi in the 
following equation 

(J c ) - 1.55((V) -(/<=)) = -3.300 log Po + Cvi , (1) 

where I c denotes the Kron-Cousins colors, (...) means 
intensity-averaged magnitudes. The above formula has been 
derived by U99 for LMC with Cvi = 15.868. Since the 
left-hand-side of Eq. (1) is reddening- free, the difference of 
CV/(LMC) and the average of this constant for the BW stars 
will yield the true distance modulus for the LMC. 

The PLC relation for the fundamental mode RR Lyrae 
stars is given by KW01 

V- 3.1(B -V) = -2.467 log P + C BV , (2) 

where the bars denote simple magnitude averages. The con- 
stant for LMC could be evaluated from the current V, B 
observations of fundamental mode RR Lyrae stars in the 
LMC by Clementini et al. (2003, hereafter C03). Because 
they give intensity-averaged magnitudes, to retain compat- 
ibility with Eq. (2), we transform these values to simple 
magnitude averages with the aid of the formulae of Kovacs 
(2002). The average of the A\ Fourier amplitude is set equal 
to 0.30, in accordance with the value obtained from the anal- 
ysis of the RRab stars in the corresponding MACHO fields 
(Alcock et al. 2003). In this way we get V = 0.032 + (V) 
and B = 0.058 + (B). With the observed average P , (V) 
and (B) of 0?58, 19.366 and 19.744, respectively, we get 
Cbv — 17.56. However, the average magnitudes computed 
by C03 for their two fields A & B differ by as much as 
0.08 mag in V, and 0.13 mag in B. They attribute this dif- 
ference to the effect of differential extinction. We checked 
this difference by computing the averages in MACHO fields 
#6 & #13, overlapping with fields A & B of C03. With a 4cr 
clipping we got (V) = 19.311 ± 0.012 (field #6, 317 stars) 
and (V) = 19.313 ±0.010 (field #13, 262 stars). With these 
small standard deviations of the means, we can state that a 
difference in the order of 0.1 mag is not probable between 
the two fields. 

In order to get another estimate on Cbv for the LMC, 
we use the data of Soszyhski et al. (2002), which cover a large 
sample of RR Lyrae stars of the Small Magellanic Cloud 
(SMC) from the OGLE database. By proceeding as above, 
from their average magnitudes and period of (V) = 19.74, 
(B) = 20.10, (Vb) = 19.45, (B ) = 19.73 and (P ) = 0.59, we 
get Cbv (SMC) = 18.01 and 17.97, depending on whether 
reddened or dereddened magnitudes are used. Although 
Eqs. (1) and (2) should, in principle be reddening-free, the 
above difference shows that the field-by-field reddening cor- 
rections applied by Soszynski et al. (2002) may not be com- 
pletely consistent with our standard extinction ratio. Since 



the method of their differential reddening correction has al- 
ready proved to be successful (at least in statistical sense, 
see U99), we adopt the value obtained by the use of their 
dereddened magnitudes. In computing the value of Cbv for 
the LMC, we have to subtract the relative distance mod- 
ulus between the LMC and SMC. From the PLC relations 
derived for Cepheids in the two clouds, U99 got a difference 
of 0.51 mag, which, with the above preliminaries, yields our 
finally adopted value of Cbv = 17.46 for the LMC. 



4 LMC DISTANCES FROM BW ABSOLUTE 
MAGNITUDES AND PLC RELATIONS 

The necessary data on RR Lyrae stars have been gathered 
from the literature during our earlier investigations (see 
KW01 and also C95, for a full list of variables subjected 
to BW analyses in earlier studies). For Cepheids we turned 
to the database of the McMaster University 2 and also to 
the radial velocity database of the Moscow University (see 
Gorynya et al. 1998). We utilized only the best data avail- 
able on fundamental mode variables. When it was necessary, 
the Cepheid I light curves published in the Johnson system 
were transformed to the Kron-Cousins system with the aid 
of the formula given by Caldwell & Coulson (1987). Occa- 
sionally, periods have been improved by using the program 
package mufran (see Kollath 1990). The variables used in 
the present study together with their derived physical pa- 
rameters are given in Table 1. (Effective temperatures have 
been computed from the cycle-averaged V — K and log<? 
values.) Three RR Lyrae stars (W Crt, SS Leo and AV Peg) 
have been excluded from the above sample in the subse- 
quent analysis, because of their discrepant positions in the 
PLC relation. For similar reasons, we also excluded X Lac 
from the Cepheid sample. 

The corresponding PLC relations are shown in Figs. 1 
and 2. The constants in Eqs. (1) and (2) have been derived 
from the BW absolute magnitudes and reddening-free col- 
ors. The straight lines represent these calibrated relations. It 
is seen that the empirical slopes of the PLC relations fit per- 
fectly the BW data (please note the highly different scales 
of the two figures). This strong correlation is very comfort- 
ing, since two, completely separate methods and datasets 
are compared (BW results, based on Galactic variables on 
one hand, and empirical results, based on LMC and globu- 
lar cluster variables on the other). The standard deviations 
of the BW points around the straight lines are 0.104 and 
0.157 for RR Lyraes and Cepheids, respectively. We recall 
that the corresponding empirical values are 0.041 and 0.058 
(see KW01 and U99). The much smaller values of these em- 
pirical standard deviations show that the BW luminosities 
still suffer from large random errors as it follows from the 
delicate nature of the BW method. 

Computation of the LMC distance modulus is straight- 
forward from the above results. It is obtained simply by the 
subtraction of the corresponding constants in the BW and 
empirical relations. By using the results of Sect. 3, from the 
RR Lyrae sample we get m — M = 17.46 + 1.09 = 18.55, 
whereas from the Cepheid sample we obtain m — M — 

2 http://dogwood. physics. mcmaster.ca/Cephcid//HomcPagc. html 
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Table 1. Derived physical parameters of the BW variables 
RR Lyraes 



Name 


logP 


T e ff 


R/Rpy 


M v 


Cbv 


SW And 


-0.354 


6702 


4.47 


0.81 


-1.09 


WY Ant 


-0.241 


6296 


6.10 


0.52 


-1.17 


X Ari 


-0.186 


6111 


6.85 


0.43 


-1.15 


RR Cet 


-0.257 


6420 


5.22 


0.75 


-0.98 


UU Cet 


-0.217 


6258 


6.25 


0.48 


-1.22 


W Crt 


-0.385 


6788 


4.90 


0.56 


-1.31 


DX Del 


-0.325 


6691 


4.62 


0.75 


-1.13 


SU Dra 


-0.180 


6293 


5.85 


0.62 


-0.94 


SW Dra 


-0.244 


6393 


5.43 


0.68 


-1.02 


RX Eri 


-0.231 


6383 


5.40 


0.70 


-0.98 


RR Gem 


-0.401 


6721 


4.60 


0.74 


-1.19 


TW Her 


-0.398 


6732 


4.55 


0.76 


-1.16 


RR Leo 


-0.344 


6458 


4.85 


0.87 


-0.93 


SS Leo 


-0.203 


6431 


6.80 


0.16 


-1.41 


TT Lyn 


-0.224 


6283 


5.72 


0.66 


-1.04 


V445 Oph 


-0.401 


6792 


4.03 


0.96 


-1.02 


AV Peg 


-0.409 


6836 


3.62 


1.18 


-0.78 


AR Per 


-0.371 


6833 


4.45 


0.72 


-1.22 


BB Pup 


-0.318 


6712 


4.40 


0.87 


-0.95 


W Tuc 


-0.192 


6316 


6.75 


0.27 


-1.26 


TU UMa 


-0.254 


6369 


5.75 


0.57 


-1.14 


UU Vir 


-0.323 


6543 


4.95 


0.75 


-1.06 


Cepheids 


Name 


logP 


T a ft 


R/R Q 


M v 


Cvi 


r\ Aql 


0.856 


5675 


60.8 


-4.04 


-2.97 


TT Aql 


1.138 


5589 


76.8 


-4.46 


-2.61 


VY Car 


1.277 


5257 


112.8 


-4.95 


-2.85 


WZ Car 


1.362 


5208 


107.6 


-4.78 


-2.38 


V Cen 


0.740 


5684 


42.8 


-3.28 


-2.60 


VW Cen 


1.177 


4998 


89.2 


-4.12 


-2.50 


XX Cen 


1.040 


5483 


64.0 


-3.96 


-2.48 


S Cep 


0.730 


5765 


44.0 


-3.42 


-2.71 


XCyg 


1.214 


5291 


100.0 


-4.72 


-2.77 


X Lac 


0.736 


5815 


55.6 


-3.98 


-3.18 


CV Mon 


0.731 


5689 


40.0 


-3.14 


-2.59 


T Mon 


1.432 


5096 


145.6 


-5.31 


-2.83 


UU Mus 


1.066 


5512 


72.0 


-4.24 


-2.65 


U Nor 


1.102 


5696 


80.4 


-4.67 


-2.89 


BF Oph 


0.609 


5903 


36.8 


-3.16 


-2.77 


BN Pup 


1.136 


5475 


74.0 


-4.26 


-2.53 


RY Sco 


1.308 


5312 


98.8 


-4.72 


-2.52 


BB Sgr 


0.822 


5645 


50.0 


-3.58 


-2.70 


U Sgr 


0.829 


5741 


51.6 


-3.75 


-2.82 


WZ Sgr 


1.339 


5064 


125.6 


-4.94 


-2.72 


RY Vel 


1.449 


5363 


132.0 


-5.41 


-2.66 


RZ Vel 


1.310 


5210 


122.8 


-5.07 


-2.88 


SW Vel 


1.370 


5231 


114.0 


-4.93 


-2.53 


T Vul 


0.647 


5986 


38.4 


-3.32 


-2.67 


U Vul 


0.903 


5964 


55.2 


-4.10 


-2.90 



15.87 + 2.69 = 18.56. Although these true distance moduli 
are magically close to each other, we have to keep in mind 
that they both have a statistical errors of « 0.03 mag. Never- 
theless, it is clear that the often cited 0.2-0.3 mag difference 
between the BW distances derived from the two groups of 
variables no longer exists. They both yield the same distance 
within a relatively small statistical uncertainty. 

It is interesting to compare the absolute magnitudes 
obtained in the present study with those derived earlier in 
other BW works. With the 18 Cepheids common with the 
sample of Gieren at al. (1998), we get and average differ- 
ence (ours-theirs) of —0.003 mag, with a standard devia- 




-.4 -.3 -.2 

Log P 

Figure 1. PLC relation for the 19 fundamental mode RR Lyrae 
stars listed in Table 1. The straight line corresponds to Eq. (2) 
with the adjusted constant given in the upper right box. 



G -5 
I 

>, 

10 -6 

in 



i-h -7 




1 1.2 

Log P 



1.4 



Figure 2. PLC relation for the 24 fundamental mode 8 Cephei 
stars listed in Table 1. The straight line corresponds to Eq. (1) 
with the adjusted constant given in the upper right box. 



tion a of 0.160 mag. A similar comparison with the absolute 
magnitudes given for p — 1.35 in Table 21 of C95, yields 
—0.162 with a = 0.094. This considerable difference results 
in a proper shift in the RR Lyrae distance modulus and 
leads to a result which is consistent with the one derived 
from Cepheids. The fact that Gieren et al. (1998) derives a 
0.09 mag lower distance for the LMC, in spite of the same 
average luminosity, is suspected to be caused by the differ- 
ence between the PLC we use and the PL they employ to 
compute the LMC distance. 

We find also a good agreement if we compare our radii 
with those of Gieren et al. (1998). For the 18 common vari- 
ables we get an average difference (ours-theirs) of —0.4 Rq 
with a = 6.3 Rq. In terms of the relative radii these figures 
become —0.04% and 6.9%, respectively. 

We can also make a comparison with the results of long- 
baseline interferometric measurements for Cepheids. For the 
two variables S Cep and r\ Aql common with the samples of 
Nordgren et al. (2002) and Lane, Greech-Eakman & Nord- 
gren (2002), we find a rather close agreement between the 
radii. By applying the necessary corrections due to the dif- 
ferent p factors, from the interferometric/radial velocity re- 
sults we obtain the radii of 42 x 1.35/1.31 = 43/?q and 
62 x 1.35/1.43 = 58-Rq for S Cep and rj Aql, respectively. 
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These values are very close to our corresponding figures (see 
Table 1). 

It is important to note that although all direct estimates 
based on double-mode variables yield LMC distance mod- 
ulus in agreement with the present result (Kovacs 2002), 
the indirect distances derived from globular cluster double- 
mode RR Lyrae (RRd) stars become higher by w 0.1 mag, if 
we use the same PLC for the LMC as in the present paper. 
Nevertheless, these revised distance moduli computed from 
the RRd populations of M15, M68 and IC4499 are still in 
the range of 18.54-18.61, which range is consistent with the 
value based on the present BW analysis. 

We mention the recent effort of Cacciari et al. (2000) 
in revisiting the BW analysis of RR Lyrae stars with up- 
dated physics. From their analysis of two stars they con- 
cluded that the derived luminosities were in agreement with 
the earlier results (meaning the continuance of the old dis- 
crepancy, however, see C03 for a different conclusion from 
the same work). 



5 CONCLUSIONS 

We have shown that by using exactly the same method and 
Baade- Wesselink parameters (e.g., velocity projection fac- 
tor), fundamental mode RR Lyrae and Cepheid variables 
yield the same distance modulus for the Large Magellanic 
Cloud. This is in variance with the earlier common belief 
that the two groups of stars lead to significantly different 
distance estimates for the LMC. The higher luminosity ob- 
tained for the RR Lyrae stars in the present study is at- 
tributed to: (a) the more accurate consideration of phase- 
dependent bolometric correction, (b) the more precise com- 
putation of the temporal gravity, (c) the better numerical 
treatment of both the grids of stellar atmosphere models 
and the observed light and velocity curves. The good agree- 
ment between the Cepheid luminosities obtained in this and 
in the earlier works suggests that the above conditions have 
already been met in those studies. 

Although clarification of important theoretical prob- 
lems (such as the accurate conversion of the projected veloc- 
ity to the photospheric value, or the application of dynamical 
atmosphere models rather than static ones) are still needed, 
and better agreement between various photometric data on 
LMC is still demanded, we think that the present result is 
encouraging and indicates the possibility of making progress 
either by increasing the sample size by more precise data, 
or by building dynamical atmosphere models. 
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